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Abstract. Analytical formulae for the calculation of secondary particle yields in p-A interactions are given. 
Fhese formulae can be of great practical importance for fast calculations of neutrino fluxes and for designing 
new neutrino beam-lines. Fhe formulae are based on a parameterization of the inclusive invariant cross 
sections for secondary particle production measured in p-Be interactions. Data collected in different energy 
ranges and kinematic regions are used. The accuracy of the fit to the data with the empirical formulae 
adopted is within the experimental uncertainties. Prescriptions to extrapolate this parameterization to 
finite targets and to targets of different materials are given. The results obtained are then used as an input 
for the simulation of neutrino beams. We show that our approach describes well the main characteristics of 
measured neutrino spectra at CERN. Thus it may be used in fast simulations aiming at the optimisation 
of the long-baseline neutrino beams at CERN and FNAL. In particular we will show our predictions for 
the CNGS beam from CERN to Gran Sasso. 
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1 Introduction 

A renewed interest in the predictions of absolute flux, 
energy distribution and flavour composition of neutrino 
beams produced at accelerators is motivated by the long- 
baseline neutrino beams in construction at CERN (CNGS 
0) and FNAL (NuMI @) and derived from the decays in 
flight of mesons generated by 400 GeV/c and 120 GeV/c 
proton interactions on graphite, respectively. 

The available experimental data on particle yields from 
high energy proton interactions on light nuclei targets, 
such as the ones used in the production of neutrino beams, 
are not extensive and one is often faced with the problem 
of extrapolating to different target materials and shapes 
and to different incident proton energies. 

Motivated by the need to estimate the neutrino flux at 
the West Area Neutrino Facility at CERN || , a consistent 
set of measurements of particle yields on Beryllium targets 
of different lengths were performed at CERN by the NA20 
Coll. f§ and subsequently by the NA56/SPY Coll. §,§, 
with 400 GeV/c and 450 GeV/c incident protons, respec- 



tively. Secondary particles (7r H 



K + , K , p, p) were 



measured in a transverse momentum range from to 600 
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McV/c, with the NA20 measurements in the high momen- 
tum range (60 GeV/c < p < 300 GeV/c) and with NA56/ 
SPY extending to the low momentum region (7 GeV/c 
< p < 135 GeV/c), thus covering also the relevant kine- 
matic range of secondary particles for the planned long- 
baseline (LBL) neutrino beams. 

In this work, in order to make these measurements of 
general applicability, the measured single-particle yields 
on Beryllium have been converted to single-particle in- 
variant cross sections, correcting for target efficiency and 
tertiary particle production. The resulting data has been 
empirically parameterized as a function of the transverse 
momentum (j>t) and the scaling variable xr = E* /Ei^ ax , 
defined as the ratio of the energy of the detected particle in 
the centre-of-momentum frame and the maximum energy 
kincmatically available to the detected particle, with a 
formula based on general physical arguments. The choice 
of the variables used in the description of the invariant 
cross section is motivated by the approximate scaling be- 
haviour of the invariant cross section, when expressed in 
terms of these variables. A satisfactory comparison of this 
empirical parameterization, based on data collected with 
400 GeV/c and 450 GeV/c protons, with invariant cross 
section measurements with 100 GeV/c protons shows that 
such approach might be useful also for the prediction of 
the NuMI neutrino beam at FNAL. 

Extension of the parameterization to target materi- 
als other than Beryllium is possible by the known depen- 



2 



M. Bonesini et al.: On particle production for high energy neutrino beams 



dence of the invariant cross section on the atomic number 
as shown in the measurements of Barton et al. 0. This 
is particularly relevant for the planned long-baseline neu- 
trino experiments at CERN and FNAL, both planning to 
use a graphite target. 

Such parameterization of the invariant cross section 
for proton-nucleus interactions could be implemented in 
a full Monte Carlo simulation of the target geometry. On 
the other hand, we show that it is possible to compute an- 
alytically particle yields from a given target by taking into 
account the geometry of the target and the contribution 
of cascade processes. Based on the NA56/SPY measure- 
ments of particle yields from targets of different lengths, 
an empirical formula describing tertiary particle produc- 
tion is derived. 

This approach, coupled to the tracking of the produced 
secondary hadrons through the neutrino beam-line ele- 
ments (magnetic lenses, collimators and decay tunnel), has 
been implemented in a fast Monte Carlo for the explicit 
calculation of neutrino beam fluxes. Effects due to reinter- 
actions in the beam-line material have been parameterized 
with good approximation. A comparison with the neutrino 
and anti-neutrino spectra measured by the CHARM II ex- 
periment at CERN is shown. Predictions for the planned 
CNGS long-baseline neutrino beam are given. 

With respect to the standard approach using complex 
Monte Carlo simulations based on hadronic cascade codes, 
such as GEISHA @ or FLUKA92 [|^s implemented in 
the GEANT packaged, CALOR95 the stand-alone 
version of FLUKA fljor MARS Q , a simple empirical 
parameterization of particle yields as described here has 
the advantage of a more transparent physical input to the 
calculations. Moreover, when coupled to the tracking of 
the secondary hadrons from the target, it becomes a fast 
neutrino Monte Carlo generator with the obvious advan- 
tage of enhanced statistics. Detailed studies of systematic 
effects due for example to horn misalignments in the neu- 
trino beam-line are made possible even with limited com- 
puter time. 



2 Parameterization of p-Be data 

An empirical formula for particle production, based on 
the yields from thick Be targets measured by NA20, has 
long ago been proposed by Malensek |Q. This has al- 
ready offered an alternative approach to neutrino beam 
computations [pi . Extensions of this parameterization to 
include effects related to the target geometry have been 
proposed [jl6|. Malensek's parameterization, however, fails 
to describe the particle production in the low momentum 
region covered by NA56/SPY measurements and cannot 
be stretched to consistently fit the whole set of available 
data on p-Be interactions. An alternative parameteriza- 
tion of NA20 data on particle yields, proposed by NA20 
gives an even less accurate description of the NA56/ 
SPY data. 

At variance with these works, the particle yields mea- 
sured by the NA20 and NA56/SPY collaborations have 



been converted to inclusive invariant cross sections for par- 
ticle production, which have then been parameterized. 



2.1 The data 

In Ref. ||, inclusive invariant cross sections for particle 
production in the forward direction in p-Be interactions 
have been reported by the NA56/SPY Collaboration. We 
have used the same method discussed in that paper to 
derive inclusive invariant cross section at all energies and 
angles from the published NA20 j| and NA56/SPY |],§ 
data on particle yields in p-Be interactions. 

The inclusive invariant particle production cross sec- 
tion is related to the measured yields from targets of finite 
length by: 



Ex ^rr = (100 -Y) 
dp J 



E 



A 



p 3 N oP X P F{L) 



(1) 



where Y is the yield per incident proton per sr (Ap/p 
%); A, p, L the atomic mass number, density and length 
of the target; Nq the Avogadro's number; E and p the 
energy and momentum of the secondary particle, respec- 
tively. The target efficiency F{L) can be estimated assum- 
ing that the produced secondary particles if absorbed in 
the target do not generate additional particles. Under this 
approximation and for forward production, one has 



F(L) 



dz 

An 



e - z /\ Pe -(L-z)/X 3 



(2) 



representing the convolution of the probability that the 
primary proton interacts between z and z+dz in the target 
and that the produced secondary particle escapes from 
the target (A p and A s are the effective mean free paths of 
primary and secondary particles respectively). This model 
neglects the production of tertiaries. This effect can be 
accounted for, in a model independent way, by estimating 
the invariant cross sections for various target lengths and 
then extrapolating the results to zero target thickness || . 

This procedure also allows to estimate at each momen- 
tum the contribution to the total yield from finite targets 
due to tertiary particles. 

For particle production at angles different from zero, 
the method is easily generalised by taking into account the 
target geometry in eq. (Q). At angles different from zero, 
however, the NA56/SPY collaboration has measured the 
production yields only at a fixed target length. In order 
to derive inclusive invariant cross sections from those data 
points, we have made the additional assumption that the 
correction for tertiary production be angle independent. 
This has been verified to be approximately true in the 
kinematic region covered by NA20 data, collected with a 
planar target setup identical to the one adopted in the 
NA56/SPY experiment. 

In this way we have derived the single-particle inclu- 
sive invariant cross sections at all available energies and 
production angles from NA20 and NA56/SPY data. Table 
[j] summarizes our results. Errors include the experimental 
errors and the uncertainty on the extrapolation procedure. 
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111.1±7.0 
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33.8i2.0 



64.5i4.7 
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10.7i0.7 
4.7i0.3 
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10.5i0.7 



2.4i0.2 
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0.300 0.300 
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64.0i4.6 
68.1i4.7 
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17.1il.2 
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14.1i0.9 
8.4i0.5 
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4.13i0.28 
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43.1i3.0 
32.0i2.1 
20.4il.7 



0.493i0.030 
0.343i0.025 
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0.114i0.008 0.152i0.010 (1.36i0.11)xl0~ 3 24.0il.5 



(1.32i0.07)xl0~ 4 
(5.57i0.52)xl0 -5 



Table 1. Inclusive invariant cross-section for pBe — » hX production as a function of Xhab — p I Pine and pr, where p and pi nc 
are the momenta of the detected particle and of the incident proton in the laboratory reference frame. Data on single-particle 
inclusive production yields of Ref. and Ref. 0], with incident proton beams of 450 GeV/c and 400 GeV/c respectively, have 
been used. 



2.2 The parameterization of inclusive invariant cross 
sections 

One of the goals of this analysis is to adopt a simple func- 
tional form for inclusive particle production which will be 
appropriate for extrapolation to different centre of mass 
energies and/or secondary particle momenta. 

Feynman has speculated in very general terms about 
the shape and energy dependence of inclusive processes 
flul , suggesting an approximate scaling behaviour of the 
single-particle inclusive invariant cross section when ex- 
pressed in terms of the transverse momentum (j>t) an d 
the longitudinal variable xf — 2p* L / \/s, where p* L and y/s 
are the longitudinal momentum of the detected particle 
and the total energy in the centre-of-momcntum frame. A 
factorization in xf and px of the invariant cross section 
has also been advocated as an experimental fact. 

An alternative scaling variable xr = E*/E^ nax , de- 
fined as the ratio of the energy of the detected particle 



in the centre-of-momentum frame and the maximum en- 
ergy kinematically available to the detected particle, was 
suggested by Yen 18 and Taylor et al. |19| and shown 
to greatly extend the range of validity of scaling at sub- 
asymptotic energies. As discussed in section 4.2, in the 
comparison of NA56/SPY and NA20 data to pBe data 
collected at lower energies, we have not observed an im- 
proved scale invariance when xr is used. Nonetheless, we 
have adopted this variable because of two practical ad- 
vantages, which translates in a simpler analytical param- 
eterization: xr is always positive by construction and can 
never be zero, except for infinite energy in the centre-of- 
momentum frame. 

Early phenomenological analyses of pp data, based 
on the constituent quark model, helped clarify that in 
hadronic processes the produced particles reflect the mo- 
tion of the constituents pp| . In particular, in the frag- 
mentation region of the projectile, the longitudinal mo- 
mentum distribution of produced hadrons reflects the mo- 
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Table 2. Values of the parameters corresponding to the best-fit of our empirical parameterization of ir^ and K inclusive 
production in p-Be interactions. Best-fit results on proton and anti-proton production data are also given (see text for details). 



mentum distribution of valence quarks inside the incident 
hadron. This suggests an x dependence of hadron produc- 
tion at large x n with a functional form similar to the one 
describing parton distributions. The functional shape of 
non-direct formation processes at small x, in which sea 
quarks are involved, however, is not easily described in 
this framework. 

After some trials and considerations of the physical 
process, the following empirical parameterization of the in- 
clusive invariant cross sections for positive sign secondary 
meson (tt + , K + ) production in p-Be interactions has been 
found to give a suitable description of data: 

(E x ^) = Ail - x R ) a (l + Bx R )x- R x 

(1 + a'(x R ) PT + b'(x R )p 2 T )e~ a '^^ (3) 

where a'(x R ) = a/x R and b'(x R ) = a 2 /2x R . 

These formula assumes an approximate factorized scal- 
ing form in x and pt- The (1 — x) a behaviour at large x 
is theoretically motivated on the basis of quark counting 
rules [^l],|2^]. The x~@ behaviour empirically accounts for 
the non-direct hadron formation mechanism at small x. 

The pt behaviour is modelled with the known expo- 
nential fall of soft interactions and a polynomial behaviour 
to interpolate the low pt part of the spectrum. The x de- 
pendence of a'(x) and b'{x) is introduced to parameterize 
the violation of pt invariance observed in the data. Mod- 
els based on the parton structure of the hadrons predict 
a Pt™ dependence of the cross section at large pt, where 
hard scattering processes take over. A reasonable param- 
eterization of p-Be data with this functional form has not 
been found. The possibility that the proposed parametric 
form fail to describe particle production in the pt region 
not covered by present data is acknowledged. 

The ratio r of positive to negative data (n + /n~ or 
K + 1 ' K~) has been instead parameterized with the empir- 
ical formulae: 

r(7r) =r - {l + x R ) r i (4) 
r(K) = r • (1 - xrY 1 (5) 

The shape of these ratios is supported by the phenomeno- 
logical analysis of pp data of Ref. p(| , showing that r(w) ~ 
1 for x ~ and rises to about 5 for x — ► 1, closely fol- 
lowing the u/d ratio of valence quarks in the projectile 

1 At large xf, xf and xr are equivalent. 



proton, while r(K) has a (1 — x)~ 3 behaviour for x — * 1. 
NA56/SPY and NA20 data only cover the fragmentation 
region of the proton at large x and the central region. At 
large x a functional behaviour similar to the one exhibited 
by pp data is expected. 

In order to keep the number of free parameters lim- 
ited, positive and negative mesons are assumed to have 
the same pr distributions. This has long been known to 
be only approximate in pp data |2^] . 

Table || summarizes the results of our best-fit to the 
data. As indicated in the table, some of the parameters 
have been fixed in the fitting procedure, since they ap- 
peared to be redundant. In the fit, 5 = 2j was chosen 
and B was set to zero, since its fitted value was found to 
be consistent with zero within errors. 

The comparison between the empirical parameteriza- 
tion and the experimental data is shown in figure |lj. The 
accuracy of the parameterization of the tt^ and if data 
is displayed in figure ^|, showing the relative discrepancy 
between our parameterization and the experimental data 
as a function of x R . The proposed parameterization gives 
an accurate fit of data with a reduced \ 2 — 0.85, while 
the reduced x 2 is somewhat larger (x 2 /Ndof — 77.1/37) 
for data. This partly reflects our difficulty to param- 
eterize the x dependence of the pt distribution, although 
about 1 /3 of the x 2 is contributed by the two data points 
(one measured by NA20 and the other by NA56/SPY) 
for positive pion production in the forward direction at 
x = 0.3, which are about 30-40% off the best-fit predic- 
tion. A reduced x 2 around 1 is obtained, if a relative error 
of 10% on each experimental point is assumed. We con- 
clude that the proposed empirical formulae are adequate 
to describe NA20 and NA56/SPY data with a 10% accu- 
racy. 

Proton and anti-proton data, which are of less direct 
interest for neutrino beams, have been parameterized with 
the following empirical formulae: 

(E x ^) pBe ^ pX = A(l + Bxr)(1 - x R ) b rt x 

(1 + ap T + YPr)e- apT (6) 

(E X -rj)p B e^pX = A(l - xh^Xx 13 X 

(l + a PT + ^p 2 T )e- a ^ (7) 
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Fig. 1. Invariant cross section as a function of pr' (top-left) positive pions; (top-right) negative pions; (bottom-left) positive 
kaons; (bottom-right) negative kaons. Data collected at the same XLab = p/ptnc, where p and Pi nc are the momenta of the 
detected particle and of the incident proton in the laboratory reference frame, are displayed with the same symbol. The best-fit 
obtained with the parameterization described in the text is superimposed. 
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Fig. 2. Percent difference between the best- fit prediction of the parameterization proposed in this work and experimental data 
(fit accuracy) for pions (left) and kaons (right) as a function of xr. 



For anti-protons a functional similar to the one given in ([j]) 
has been adopted, except that an exact factorization in xr 
and pt has been assumed, since this was sufficient to give 
a reasonable fit to data. For protons the "leading particle 
effect" had to be taken into account. A reasonable fit to 
data has been obtained by following the empirical observa- 
tion that the longitudinal momentum distributions of the 
leading nucleon in pp collisions is flat pi| , which translates 
into a linear rise of the inclusive invariant cross section as 
a function of xr. The transverse momentum distribution 
is also affected by the leading particle effect, resulting in 
an enhanced leading particle production in the forward 
direction (see for example |l^,^3|). In the proton fit, this 
is empirically accounted for by the term (I — xji) hpT . The 
parameter a, that controls the shape of the px distribution 
for non-leading particle production, has been assumed to 
be the same for protons and anti-protons. 

Results of these fits are also given in Table [2] and the 
comparison between the empirical parameterization and 
the experimental data is shown in figure ||. Our param- 
eterization gives a satisfactory description of proton and 
anti-proton inclusive production, with a reduced \ 2 about 
1, in the range covered by NA56/SPY and NA20 data. At 
larger values of x, proton production by means of diffrac- 
tion should occur. This is not described by our parame- 
terization. 



3 Production of neutral Kaons 

The knowledge of neutral kaon production in p-Be inter- 
actions is important for the exact calculation of the v e 
background in neutrino beams coming from K® L decays. 
As an example, in the case of the WANF beam serving 



the CHORUS and NOMAD experiments, the estimated 
contribution to the v e content of the beam due to K\ 
is around 15%. In addition, the knowledge of the neutral 
kaon production is also of interest for the experiments on 
CP-violation using neutral kaon beams and it is particu- 
larly relevant in searches for rare neutral kaon decays. 

A rather complete set of data on Kg production at 
angles different from zero was collected by Skubic et al. 
p6[ , using a 300 GeV/c proton beam hitting Be targets. 
Measurements of K s production in the forward direction 
from 200 GeV/c protons on beryllium were performed 
by Edwards et al. p7y . Both these experiments, however, 
covered a region of secondary momenta corresponding to 
xf > 0.25, which only partially overlaps to the momentum 
region important for present and future neutrino beams. 

Additional information on neutral kaon production can 
be obtained starting from the available measurements of 
charged kaon production in p-Be interactions. Charged 
and neutral kaon production rates can be related in the 
frame of a simple quark parton model, in which valence 
quarks (q v ) and sea quarks (q s ) are considered. 

In the simplest model, isospin symmetry may be as- 
sumed, giving u v /d v = 2, u s — u s = d s — d s and s s = s s , 
from which one finds that the production of charged and 
neutral kaons should be related by 

N(K° S ) = N(K° L ) = \{N K+ + 3 ■ N K -) (8) 

This model has been shown to provide an accurate de- 
scription of the asymmetry of K° / K° production in the 
range 0.18 < xp < 0.36 in neutral kaons beam-lines p5[ . 
We have verified that it also agrees within 15% with direct 
measurements of K s production up to around xf = 0.5. 
At larger values of xf , i.e. in the fragmentation region of 
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Fig. 3. Invariant cross section as a function of pr- (a) protons; (b) anti-protons. Data collected at the same x^at = p/pi nc , 
where p and pine are the momenta of the detected particle and of the incident proton in the laboratory reference frame are 
displayed with the same symbol. The empirical parameterization described in the text is superimposed. 



the proton, we find that a reasonable description of ex- 
isting K s production data requires a "dynamic" model in 
which the u v /d v ratio is let to be ^-dependent and es- 
timated from the ir + /tt~ production ratio, following the 
arguments given by Ochs po| , who has demonstrated a 
remarkable empirical similarity between the n + pro- 
duction ratio in pp collisions and the u/d quark ratio mea- 
sured in deep-inelastic lepton-proton scattering. The ex- 
ploitation of this model is beyond the aim of this paperf] 
and neutral kaon production is always estimated from eq. 
@ in the following. 



4 Particle production for neutrino beams 

4.1 Scaling to targets of different materials 

Beryllium targets have been used in most neutrino beams 
derived from proton beams extracted from the Super Pro- 
ton Synchrotron (SPS) at CERN, but the use of graphite 

2 A similar "dynamic" prediction of neutral kaon production 
was already considered in ref . . These models however seems 
less accurate than formula Q in the description of the K° /K° 
production asymmetry at small/intermediate xf. 



as target material seems promising in view of operation 
with short spills in fast extracted proton beams. In partic- 
ular, the use of a graphite target is foreseen in the design of 
both the CNGS and the NuMI beams Prescriptions 
to rescale the inclusive invariant cross-sections to different 
target materials arc given in the following. 

Invariant cross sections E d ^ p3 for hadron- nucleus in- 
teractions (hA i— » h'X) depend on the mass number A of 
the target nucleus, via parameterizations of the type: 




where a value for a — 2/3 would correspond to the case 
where particle production off a nucleus is identical to the 
production off a single nucleon per inelastic collision. In 
accordance with the scaling hypothesis, a has been found 
to be weakly dependent on the incident beam momen- 
tum. It depends on the incident hadron type h and it is a 
smooth function of px and if of the produced hadron. It 
has moreover been experimentally observed that, in first 
approximation, it is independent of the detected particle 
type, with perhaps the exception of anti-protons H. 

A parameterization of a as a function of Xf has been 
proposed by Barton et al. on the basis of their and pre- 
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vious data at px — 0.3 GeV/c M. A px dependence of 
a has been clearly observed in Kg and A production by 
Skubic et al. |26| . A suitable representation of the whole 
set of data can be obtained with the parameterization: 

a(x F ) = (0.74 - 0.55-a; i r+0.26-a;|)-(0.98+0.21-p|) (10) 

where the x f dependence is taken from the fit of Barton 
et al. at px = 0.3 GeV/c and the px dependence is fitted 
to Skubic et al. data and normalized in such a way that it 
reduces to the parameterization of Barton et al. at px = 
0.3 GeV/c. 

In conclusion, a conservative estimate of the uncer- 
tainty in the extrapolation from beryllium to carbon data 
sits (in the px range of interest for neutrino beams: up 
to ~ 600 MeV/c) around 5%, on the top of a measure- 
ment error of 5-10%, depending on the secondary momen- 
tum, for the cross sections on beryllium. The estimate of 
this systematic uncertainty is based both on data collected 
by Barton et al. and on the extensive compilation of 
J. Kuhn on nuclear dependence for pA — > n~ X interac- 
tions M. 



4.2 Scaling to different centre of mass energies 

The NuMI neutrino beam at FNAL is planned to be de- 
rived from a primary proton beam of 120 GeV/c momen- 
tum, resulting in a centre-of-mass energy about two times 
smaller than that available at NA56/SPY and NA20. Be- 
sides its general interest, a test of the scaling hypothesis of 
one-particle inclusive invariant cross-sections is thus rel- 
evant to assess the range of validity of the proposed pa- 
rameterization. This has been studied by comparing our 
parameterisation to available pA data collected at differ- 
ent centre-of-mass energies. 

In figure |^, our prediction for one-particle inclusive in- 
variant cross section are compared to data collected with 
100 GeV/c protons on carbon target JtJ, where the ex- 
trapolation from Be to C has been made using formulae 
(^) and (|To|). As discussed above xr — E*/E^ ax has been 
used as scaling variable. The agreement between pion data 
shown in the figure and our empirical parameterization is 
excellent up to about x ~ 0.8. A good agreement is also 
found when kaon data are considered, although the preci- 
sion of kaon data from ref. R is poorer. 

Figure [| shows a comparison of our parameterization 
with one-particle inclusive production data measured with 
24 GeV/c protons hitting a Beryllium target |29]. In that 
work, Lorentz invariant particle densities oj(j>Lab, 9) were 
measured as a function of the particle momentum pxab 
and production angle 9. These data have been converted 
into invariant cross sections using the relation 



d 3 a 1 
E x — - = -a abs uj[p Lab ,9) 
dp 6 2 



(11) 



distributions for 7r + and K + , although the estimated pro- 
duction of 7r + is about 35 ± 15% lower than that mea- 
sured. This is also true for negative pions, not displayed 
in the figure, while the agreement is somewhat worse for 
the other particles. 

Given its general interest, we have tested the scaling 
hypothesis on the same data also using xp and xxab = 
p/Pinc, the latter defined as the ratio of the momentum 
of the detected particle in the laboratory reference frame 
to the momentum of the incoming proton, as longitudi- 
nal variables. A better agreement with 7r ± data is ob- 
tained, if XLab is used, which however shows a worse agree- 
ment to kaon data^j. The variable Ay — y max — y, where 
y = \ log((-E + px)/{E — px)) is the rapidity of the pro- 
duced particle and y m ax is the maximum rapidity kine- 
matically available to that particle, has also been tried. 
At variance with the previuos longitudinal variables, Ay 
is Lorentz invariant and scales the phase space at differ- 
ent centre-of-mass energies independently of the reference 
frame. Although we find a somewhat improved agreement 
to data in this case, scaling to data collected with 24 
GeV/c protons is only approximate even with this vari- 
able. 

Still, in the centre-of-mass energy range of interest for 
present high-energy neutrino beams, the agreement be- 
tween our parameterization scaled according to xr and 
the data is satisfactory. 



4.3 Yields from finite length targets 

With reference to eq. (0), the differential particle produc- 
tion along the target can be parameterized as: 



dY{E, PT ,z) _ N pX p p 3 



dz 



100 A E 



(12) 



where f(z)dz is the probability that the outgoing particle 
be produced at a depth z to z + dz inside the target and 



the other quantities have been introduced in Section 2.1 



In general f(z) will also depend on the production an- 
gle of the secondary particle |lq|. In the naive reabsorption 
model introduced in Section |2.lL one has: 



f(z,9) = -U 

Ap 



-z/\ Pe -z'{z,e)/\, 



(13) 



which can be obtained by differentiating the target effi- 
ciency F(L) of equation (||), with z'(z, 9) representing the 
target length that the produced hadron has to cross to es- 
cape from the target [L — z for forward production) . This 
naive model is only adequate to give a first order descrip- 
tion of particle production from relatively thin targets or 
at large x, but in general particle production by means of 
cascade processes cannot be neglected. 

Since all types of hadrons can be produced when an 
energetic hadron of any flavour interacts in the target, the 



where a a bs is the absorption cross section for pBe fl31fl . 
A reasonable agreement is observed in the shape of the 



3 All the variables are off roughly by the same amount in the 
description of the K/n production ratio. 
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production of particles from a target of finite length is gov- 
erned by a set of coupled transport equations that depend 
on the properties of the particles, their interactions and on 
the structure (geometry and material) of the target. Given 
their complexity, the solution of these equations is in gen- 
eral addressed numerically with hadron cascade/transport 
codes, calibrated on a large set of hadron production data. 
However, in the spirit of this work, the main features of the 
cascade process in needle-shaped targets, like those used 
for neutrino beams, can be described in a simple para- 
metric form and tuned on NA56/SPY data on yields from 
targets of different length and geometry. 

In hadronic interactions, about half of the available en- 
ergy is typically dissipated in multiple particle production 
(inelasticity k ~ 0.5) with average transverse momentum 
(pt) — 0.35 GeV/c . The remainder of the energy is car- 
ried by fast forward-going "leading" particles. Only these 
particles are responsible for the propagation of the cascade 
in thin and long targets, since they are almost collinear 
with the primary beamF] and their energy is large enough 
to result in a sizeable yield of additional energetic parti- 
cles. 

In the collinear approximation and assuming that lead- 
ing particles in pA interactions are mostly protons, thus 
characterized by an effective mean free path A p equal to 
that of primary protons, the naive reabsorption model of 
equation (fL3|) can be improved by the expansion: 



(14) 



f(z) = + A h {x)^-]e-*'W x 

A 



A, 



where the second term in brackets accounts for particle 
production by reinteractions of secondary particles (ter- 
tiary particle production) , while higher rank contributions 
are neglected, since neutrino targets are typically of order 
2A P and the less energetic is the reinteracting particle the 
lower is the yield of produced particles. 

In equation (|l4]), Ah(x) weights the probability that 
reinteractions of secondary particles will result in a pro- 
duced hadron h of fractional longitudinal momentum x 
and is given by: 



A h (x) 



h' 



dx'G h ,^ h (x,x')$ h .(x') (15) 



where Gh>—,h{x,x') is the px integrated cross section for 
inclusive production of h in h' interactions and <2v(x') is 
the flux of h! particles. Thus, Ah(x) is a function of x and 
h only. Its dependence on these variables has been derived 
from NA56/SPY measurements of inclusive yields in the 
forward direction from finite length targets. 

In the short target approximation, valid until the tar- 
get length is L << XpX s / (A s — A,,), integration of equation 
( |l4| ) in the forward direction (i.e. z' — L — z) predicts the 
fraction t(L) of tertiary particle production to be approx- 



4 The divergence of leading particles with respect to the pri- 
mary beam is typically (p T )/(l — k)p inc ~ few mrad for proton 
energies of order a few hundred GeV. 
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Fig. 6. Fraction of tertiary particle production as a function 
of the target length as derived from published NA56/SPY data 



imately a linear function of L given by: 



t(L) ~ A h (x) 



L 

2A„ 



(16) 



This linear behaviour is in agreement with experimental 
data, as shown figure ^, where the fraction of tertiary pro- 
duction as a function of target length for 7r ± and .fT as 
derived from NA5 6/S PY data is displayed, and was ex- 
ploited in section 2.1 to derive inclusive invariant cross 
sections. 

Compared to the naive absorption model, the excess 
of particle production present in the NA56/SPY data in- 
creases with decreasing momentum. Figure shows the 
fraction of tertiary particles produced in the forward direc- 
tion for a 100 mm target as a function of the fractional lon- 
gitudinal momentum XL ao — p Lab I Vine- The curves shown 
in the figure represent the best-fit of Ah(x) to data ac- 
cording to the empirical parameterization: 



A h (x Lab ) = A{ er (l-x Lab ) b 



(17) 



where a common value of A 1 ^ and &^ er has been chosen 
for positive and negative particles. The best-fit results are 
reported in table [3| 

In this simple model, no pt dependence of tertiary par- 
ticle production is predicted. Indeed, assuming Feynman 
scaling and the (approximate) factorization in x and px of 
the invariant cross-section, the pt distribution of tertiary 
particles is expected to be the same as that of secondary 
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Fig. 7. Fraction of tertiary particle production in a 100 mm 
Be target as a function of XLab- Data are derived form pub- 
lished NA56/SPY data [0] and the curves represent the best-fit 
described in the text. 



A ter 


Ater 


b ter 
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°K 


0.80 


1.56 


7.3 


10.1 



Table 3. Best-fit values for parameters of equation (|17]) to 
NA56/SPY data (see text for details). 



particles. It follows that the divergence of leading particles 
with respect to the primary beam direction is in general 
negligible as compared to the typical production angle of 
lower momentum tertiary particles ( collinear approxima- 
tion). These considerations are supported by NA56/SPY 
data that show no pt dependence of the transverse mo- 
mentum distribution of tertiary particle production in a 
300 mm Be needle-shaped target once mere geometric ef- 
fects are considered (see figure 20 of Ref. [|). 



Equation (gj) is valid for NA56/SPY and NA20 tar- 
gets, but fails to give the correct prediction for tertiary 
particle production when the target length is of order 
Xgeom = R- Pi I {pt) , where R is the transverse dimension 
(radius) of the target and pi = (1 — k) ■ Pi nc is the typical 
momentum of the leading particle. In this limit, there is a 
sizeable probability that some of the leading particles also 
escape from the side of the target before interacting. This 



can be accounted for, modifying equation (|wj) into 

dy 



m = 



l + A h {x) / yw z (y) 



o A E 



z'(z,0)/X B 



(18) 

where w z {y) is an acceptance factor, which in general de- 
pends on the pt and energy distributions of the leading 
particles and on the target geometry^. For cylindrical sym- 
metry, one can write: 



w z {y) 



(Pt) 



Pt 



h(p T )p T dp T 



(19) 



where pip ax ~ piR/(z — y) in the collinear approximation 
and h(pr) is the pr distribution of the produced particles. 

In the simulation of neutrino beams, discussed in the 
following, the computation of this acceptance factor has 
been performed under the approximation that leading par- 
ticles carry on average half of the primary beam momen- 
tum and that their transverse momentum distribution is 
represented by a pure exponential of average transverse 
momentum (pr) - For cylindrical symmetry this allows the 
analytical integration of equations (jll^) and (|l^) . In case 
of continuous targets of radius R, the integration yields: 



/(*) 



Ar, 



1 + A h (x) — (1-e 

An 



-z'(z,e)/x„ 



_ , (20) 
which reduces to equation (J14|), predicting a linear increase 
of the differential production of tertiary particles along the 
target, until z << X geom (short target approximation). At 
larger depth inside the target, the differential production 
of tertiary particles will tend to saturate to a value pro- 
portional tO Xgeom /A p . 

In case of a segmented target, made of several rods in- 
terleaved by air gaps, as that used in the WANF beam and 
considered for the forthcoming CNGS beam at CERN, the 
details of the target geometry have to be considered. The 
analytical integration of equations (|l8|) and (|l9|) is still 
possible under the same approximation as that adopted 
to derive equation (20) and it has been used in the evalu- 
ation of the efficiency of segmented targets. 

In Figure || we present predictions for secondary meson 
production at target level, based on our parameterization 
of the inclusive invariant cross-sections corrected for target 
efficiency. Momentum spectra (top) , angular distributions 
(middle) and exit point profiles (bottom) of the secondary 
mesons in the case of the CNGS target geometry and for 
an angular acceptance of 30 mrad (much larger than the 
horn aperture of the CNGS beam-line described in the 
next Section) are shown. For sake of comparison, the pre- 
dictions of a FLUKA based full Monte Carlo program | fl2"[ 
[jj are also shown. The agreement is remarkable, provided 
that tertiary particle production is taken into account in 
the evaluation of the target efficiency, as discussed in this 
Section. 



5 In equation ([18]), the collinear approximation has been 
again advocated to factorize the acceptance factor and the rein- 
teraction probability of leading particles. 
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Fig. 8. Longitudinal (upper panels) and transverse (middle panels) momentum distributions of secondary mesons from the 
CNGS target (within 30 mrad acceptance). On the lower panels the exit points of the secondary mesons from the CNGS target 
(within 30 mrad acceptance) are shown. The dotted line is the FLUKA prediction, the dashed line is our computation without 
tertiary contribution, the solid line is our prediction with tertiary mesons. 
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A slight disagreement in the px distributions of all 
charged particles is visible mainly at large pt, that trans- 
lates in a slightly more forward peaked distribution of the 
exit points along the target. This disagreement appears in 
a kinematic region not covered by NA56/SPY and NA20 
data. However, as remarked in the next Section, neutrino 
beam predictions are only slightly affected by this large pr 
region, since it falls outside the momentum and angular 
acceptance selected by the focusing optics. The systematic 
difference between the two models is thus negligible. 



5 Simulation of neutrino beams 

As it is well known, a "classic" wide band neutrino beam 
is produced from the decay of mesons, mostly 7r's and if's. 
Mesons are created by the interaction of a proton beam 
into a needle shaped target, they are sign-selected and 
focused in the forward direction by two large acceptance 
magnetic coaxial lenses, conventionally called at CERN 
horn and reflector, and finally they are let to decay into 
an evacuated tunnel pointing toward the detector position. 

In case of positive charge selection, the beam content 
is mostly from the decay of ir + and K + . Small con- 
taminations of (from the defocused 7r~ and K~) and 
v e (from three-body decay of K's and /z's) are present at 
the level of few percent. 

As an example, a schematic layout of the future CERN 
to Gran Sasso neutrino beam, CNGS Q , is shown in Fig- 
ure | Its main characteristics, relevant for beam simula- 
tion purposes, are listed in Table |] together with those of 
the old West Area Neutrino Facility (WANF) at CERN in 
the configuration set up for the CHARM II experiment ||. 

The neutrino fluxes for such a kind of beam are rela- 
tively easy to predict once the secondary meson spectra 
are known, because the meson decay kinematics is well 
understood and the geometry of the decay tunnel is quite 
simple. 

Uncertainty in the estimation of the neutrino fluxes 
could arise because secondary mesons are selected over a 
wide momentum range and over a wide angular acceptance 
(~ 20 mrad). 

Re-interactions of secondary mesons in the target and 
downstream material contribute to reduce the neutrino 
fluxes and increase the uncertainty in the calculations 
(mainly for the wrong sign and wrong flavour contamina- 
tions) . They are generally minimized using a target made 
of a number thin rods of low Z material interleaved with 
empty spaces (to let the secondary mesons exit the target 
without traversing too much material). In addition the 
amount of material downstream of the target (i.e. horn 
and reflector conductor thickness) is kept to the minimum. 

The parameterization of the secondary meson produc- 
tion from protons onto a thin target, proposed in this pa- 
per, is thus well suited to be used in neutrino beam simu- 
lations both because it extends its prediction over a wide 
range of longitudinal and transverse momenta and also 
because the small fraction of tertiary production from re- 
interactions in the target and downstream material can be 



acco unt ed for with the approximations described in Sec- 
tion o 



A comparison of the neutrino flux prediction based 
on the proposed parameterization with some measured 
spectra is thus an effective estimator of the quality of 
the secondary mesons parameterization. For this purpose, 
our parameterization has been coupled with a neutrino 
beam simulation program able to provide rapid and accu- 
rate predictions of neutrino spectra at any distance (i.e. 
short and long base line). The comparison has been per- 
formed both with already published data (CHARM II) 
and with predictions for the future CNGS long-baseline 
neutrino beam generated with GEANT and/or FLUKA 
based Monte Carlo programs. 



5.1 The simulation program 

The program is a stand-alone code developed as a tool that 
allows to vary and optimize all elements and the geometry 
(in 3-D) of the beam line providing the results in terms of 
neutrino spectra and distributions at large distance with 
high statistics and in short time. As a by-product, the pro- 
gram provides also spectra and distributions of secondary 
hadrons and muons along the beam-line and in the muon 
monitor pits after the hadron stop at the end of the decay 
tunnel. 

The underlying idea is that in order to produce rapidly 
a neutrino spectrum at large distance over a small solid 
angle (typically dfl ~ 10~ 10 for the future LBL beams 
of CERN and FNAL), one has to force all the mesons to 
decay emitting a neutrino, and force all neutrinos to cross 
the detector volume. A weight is then assigned to each 
neutrino, proportional to the probability that this process 
actually happened. 

In practice this method is implemented by subdividing 
the simulation into four subsequent steps (as described 
in detail in the following). The weight assigned to each 
neutrino is the product of factors originated in each step 
times the solid angle subtended by the detector: 



det • 



(21) 



5.1.1 Mesons production along target 

First the total number of secondary mesons (tt + , tt~ , K + , 
K~ and K®) to be generated is calculated by integra- 
tion of the yield calculated in our model over the interest- 
ing range of longitudinal and transverse momenta. These 
numbers obviously depend on target material, density and 
length and on the number of protons on target. 

Secondary mesons are generated along the target ac- 
cording to the distribution of proton interaction points. 
The latter depends on the proton beam size and diver- 
gence, on the target thickness, T tar g, and on the proton 
interaction length, A p , of the target material. Momen- 
tum and angular distributions of the mesons are sampled 
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Fig. 9. Schematic layout of the future CNGS neutrino beam line. 



Beam line 


WANF 


CNGS 


Target material 


beryllium 


graphite 


Target rod length 


10 cm 


10 cm 


Target rod diameter 


3 mm 


4 mm 


Number of rods 


10 


8 


Rod separation 


9 cm 


9 cm 


Additional end-rod length 




50 cm 


Proton energy 


450 GeV 


400 GeV 


Proton beam focal point 


50 cm from start of target 


Expected pot/year: 


< 2 x 10 19 


4.5 x 10 19 


in shared SPS mode 


in dedicated SPS mode 




7.6 x 10 19 


Horn & Reflector focusing momenta 70-70 GeV 


35-50 GeV 


Horn & Reflector length 


6.56 m 


6.65 m 


Horn & Reflector current 


100-120 kA 


150-180 kA 


Horn distance from focal point 


11.5 m 


2.7 m 


Ren. distance from focal point 


82.8 m 


43.4 m 


Horn acceptance 


~ 10 mrad 


~ 20 mrad 


Decay tunnel length 


285 m 


992 m 


Decay tunnel radius 


0.6 m 


1.22 m 


Tunnel vertical slope 


+42 mrad 


-50 mrad 


Pressure in decay tunnel 


1 Torr 


1 Torr 



Table 4. Main parameter list of the WANF (CHARM II set-up) and CNGS neutrino beam lines. 



from the proposed parameterization. The weight associ- 
ated with this step is 

Wx = {l- e - T *-' 9 /V) e -2/^ (22) 



where Ah (x) , w z (y) and the approximation related to this 
approach have been discussed above. In case of the seg- 
mented target considered for the CNGS beam, this model 
gives the results of figure ||. 



because protons are forced to interact in the target. 

Meson trajectories in the target are calculated and 
their length, z'(z, 9) used to estimate the probability that 
the mesons exit the target without re-interacting: e~ z / As , 
where A s is the meson interaction length in the target. 

Tertiary contribution is added following the parame- 
terization described in Section 4.2 : 



W 2 



-z'/\ 3 



dy 



l + A h (x) I ^w z (y) 
o A P 



(23) 



5.1.2 Meson tracking in the neutrino beam-line 

The trajectory of each meson in the beam-line is calcu- 
lated, taking into account the tracking in the magnetic 
field of horn and reflector, until it hits the walls of the 
decay tunnel or the collimators. The amount of material 
crossed by the particle is also recorded. 

The meson is forced to decay along its trajectory, traj, 
accordingly to its decay length, Xdec- The weight associ- 
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ated to this process is 

W 3 = (1 - e -* raj/Ad ")e" z - t/As 



(24) 



where Zi nt is the amount of material crossed up to the 
decay point and A s is the interaction length in that mate- 
rial. 

Contributions due to reinteractions of secondary par- 
ticles in the material along the beam line are taken into 
account with a parameterization similar to that used for 
the target but in the "short target" approximation. The 
corresponding weight is: 



W A = l + A% r (l-x)"» 



2A„ 



(25) 



as derived in formulae ( |l6| ) and ( |l7| ) of Section 4.3 



5.1.3 Neutrino production from mesons 

For each meson a neutrino is produced; its flavour and its 
momentum distribution in the parent meson rest frame 
depend on the decay mode and branching ratio, B.R.. The 
neutrino direction in the laboratory frame is determined 
requiring that it crosses the detector volume. The angle, 
9 SV , between parent meson and neutrino directions allows 
to calculate the Lorenz boost of the neutrino from the 
meson rest frame to the laboratory frame. This in turns 
allows to obtain the neutrino momentum in the laboratory 
frame. 

The weight associated to this process is proportional to 
the probability that the neutrino is emitted in the detector 
direction. This is obtained by simply boosting back the 
solid angle subtended by the detector in the meson rest 
frame, where the neutrino is emitted isotropically; 



W 5 = B.R. x (m s /(E s - P s cosfl sl/ )) 2 



(26) 



detector position and the neutrino source as it is the case 
for classical unweighted methods. 

In the classical unweighted case, only a fraction of the 
pions (typically 5-10%) decay before interacting (either 
in the beam-line material or in the decay tunnel walls); 
in addition neutrinos are spread over a wide solid angle 
(about 1 mrad in the CNGS case) because of the decay 
kinematics. It follows that, to enhance statistics, at large 
distance the neutrino spectra need to be computed on a 
surface much wider than the actual detector area, relying 
on the fact that the spectra shapes varies slowly with the 
radius. In the CNGS case an accuracy of a few percent 
can be achieved with several millions protons on target if 
a detector area larger than ~ 10 4 m 2 is considered. 

In the parameterization case, since all mesons - within 
the focusing optics acceptance - are exploited to produce 
neutrinos in the detector, the statistical accuracy is inde- 
pendent from the detector distance and proportional to 
the inverse of the square root of the number of generated 
positive pions (for beams), namely about the number 
of generated proton interactions on target. An accuracy 
better than a percent is thus obtained with less than 10 5 
p.o.t., for any size of the detector surface. 



5.3 Simulations of past and future neutrino beams 

In order to give a more quantitative appreciation of the 
accuracy that one can obtain in the simulation of neu- 
trino beams using the generator described in the previous 
Sections, we present the comparison with the published 
neutrino spectra measured with the CHARM II detector 
[|§ exposed at the CERN-WANF beam. 

In addition we present the comparison between the 
CNGS simulation based on our method and that based 
on the FLUKA stand-alone program for secondary meson 
production interfaced with GEANT for particle tracking. 



where m s , E s and P s are the mass, energy and momentum 
of the secondary meson. 



5.1.4 Neutrino production from muons 

Muons are produced in the decay of secondary mesons 
taking into account the correct kinematics (branching ra- 
tio and polarisation) . Muons are also tracked through the 
neutrino beam line and forced to decay to produce neutri- 
nos in the detector direction. An additional weight, 
(equivalent to W3 for meson decays) , is introduced for the 
neutrinos from muon decay. A weight Ws^ replaces W5 
accounting for the the muon decay kinematics (including 
polarisation) . 



5.2 Statistical accuracy 

The statistical accuracy of this way of simulating neutrino 
beams does not depend much on the distance between the 



5.3.1 Comparison with CHARM II data 

The WANF neutrino beam line at CERN, to which the 
CHARM II detector was exposed, is well described else- 
where Jp|. Its main characteristics has been summarized in 
Table [ifThe facility was run, during several years of oper- 
ation, either selecting positive charged particles (z/ M beam) 
or negative ones (V^ beam). Neutrino/anti-neutrino inter- 
actions were collected in the CHARM II detector and fully 
reconstructed [^3| . 

In figure Figure 10 -top, we show the comparison be- 



tween the measured neutrino fluxes in case of positive 
mesons focusing, [y^ beam with v ^ contamination) and 
the simulation performed with our method and based on 
10 5 p.o.t.. On the left hand side a logarithmic scale is used 
to make evident the spectral behaviour at high energy as 
well as the wrong sign contamination; on the right hand 
side a linear scale is used for a better appreciation of the 
focusing/defocusing effect. In figure Figure |To[-6o^^om, the 
same plots are shown in the case of negative meson focus- 
ing (V/j beam with contaminations). 
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Fig. 10. The WANF neutrino (top) and anti-neutrino (bottom) fluxes at the CHARM II detector: the dotted lines are experi- 
mental data from Ref. [33[ |, the continuous line is the beam simulation. On the left, logarithmic scale is used to make evident 
the spectral behaviour at high energy as well the wrong sign contamination; on the right, linear scale is used for a better 
appreciation of the focusing/defocusing effect. 
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The overall agreement is at the percent level, with at 
most 10% disagreement on a bin per bin basis. 

The high energy tails of the distributions are domi- 
nated by the production of high energy secondary mesons 
peaked in the forward direction, and are practically in- 
sensitive to the magnetic focusing. In fact high energy 
mesons, with small angular aperture, travel most likely 
through the neck of the horn where they are hardly de- 
viated. The good agreement between simulation and data 
indicates that high Xf production on target is well sim- 
ulated and that re-interactions on the material along the 
beam-line is correctly taken into account. 

In the focusing/defocusing energy range the agreement 
is an indication that low xf production is correctly gen- 
erated at least up to ~ 10 mrad (the WANF optics accep- 
tance). The fact that also the wrong sign contamination 
in the simulation behaves as the data, means that tertiary 
production in target and down-stream material (mainly 
the horn neck) is described to a sufficient level of approx- 
imation. 



5.3.2 The CNGS beam-line 

As mentioned earlier, our neutrino beam generator was 
originally developed to allow a rapid optimization of LBL 
neutrino beams. To test the reliability of it, we have per- 
formed a detailed simulation of the CNGS LBL neutrino 
beam for an extensive comparison with the full beam sim- 
ulation based on the FLUKA stand-alone package |l2| for 
secondary particle production and on GEANT for track- 
ing along the beam line [EJ . 

Apart from the tunnel geometry and the focusing op- 
tics, the main differences of the CNGS beam-line with re- 
spect to the WANF are the target material {carbon instead 
of beryllium) and the proton beam energy (400 GeV in- 
stead of 450 GeV) . These differences have been taken into 
account in our simulation with the scaling laws proposed 
and described in Section 4.2 and 4.3. 

As remarked in Section 4.3, the agreement at the target 
level between our calculation and the FLUKA based full 
Monte Carlo program is noticeable, provided that tertiary 
production is taken into account in the former case. The 
slight disagreement visible in the pt distributions of all 
charged particles mainly at large pt translates in a slightly 
more forward peaked distribution of the exit points along 
the target. 

Remarkably enough, as far as the neutrino spectra at 
large distance are concerned, the discrepancies between 
the two models do not propagate with the same strength 
(see Figure [ll]). This is because in the momentum range 
and angular acceptance selected by the focusing optics 
both particle production models reproduce very well the 
available experimental data. 

Before concluding, it is worth mentioning that early 
simulations of the CNGS beam line based on GEANT 
stand-alone disagreed with those presented here, being 
too optimistic by more than ~ 20%. On the other hand 
FLUKA stand-alone gives results fully compatible with 
those presented here. 



6 Conclusions 

Empirical formulae for single-particle inclusive invariant 
cross sections in p-Be interactions have been derived, on 
the basis of single-particle inclusive production data col- 
lected by the NA20 g] and NA56/SPY [§§ experiments. 
These formulae reproduce the experimental data within a 
10% accuracy. 

The hypothesis of Feynman scaling has been verified 
to hold with our parametrisation giving a suitable repre- 
sentation of production data collected over a wide range of 
primary proton beam energies (from 24 GeV to 450 GeV). 

Prescriptions to extrapolate this parameterization to 
finite targets and to targets of different materials have 
been given. 

The results obtained have been used as an input for the 
simulation of neutrino beams. A comparison to data col- 
lected by the CHARM-II neutrino experiment at CERN 
has demonstrated the capability of this approach to pre- 
dict the past. 

These formulae can be of great practical importance 
for fast calculations of neutrino fluxes and for designing 
new neutrino beam-lines. Thus they may be used in fast 
simulations aiming at the optimisation of the long-baseline 
neutrino beams at CERN and FNAL. Predictions for the 
neutrino spectra of the CNGS beam from CERN to Gran 
Sasso have also been given. 
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